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A number of calculations have appeared for the Bose-Einstein condensation temperature (T c ) of 
a weakly repulsive dilute Bose gas. After a short survey of previous work, I point out several issues 
related to these approaches, and outline future calculations. 
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I. INTRODUCTION 

During the recent upsurge of activity in the field of 
Bose-Einstein condensation (BEC), some older questions 
concerning the weakly interacting Bose gas have been re- 
vived. One such un-resolved question is the effect of a 
weak repulsion on the critical temperature T c . Calcu- 
lations on this question has produced widely dissimilar 
results: increases and decreases of T c proportional to a, 
y/a, a In a, etc. have been reported by various authors. 
(Here a is the scattering length.) It is only during the 
past 2-3 years that the community has been converging 
toward a consensus. 

Considering the diverse approaches to the T c problem, 
and the variety of issues that have arisen in this con- 
nection, it seems appropriate at this point to attempt to 
draw some connections between the different approaches, 
and point out some misconceptions and fine points. The 
purpose in this article is to review and comment on ex- 
isting work, rather than to present new calculations. 

The model is described by the Hamiltonian 
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E 



q-kMq 
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where ey = k 2 /2m is the free-gas spectrum, and b, w 
are bosonic operators. In principle, one is interested 
in any possible interaction function f/(k), but the sim- 
plest possible form is a delta function in real space, so 
that E/(k) = U is momentum-independent. To first or- 
der U is related to the s-wave scattering length a by 
U = 4:Trh 2 a/m. A dimensionless measure of the interac- 
tion is the quantity art 1 / 3 , where n — N/V is the density. 
The program is to examine the shift in the critical tem- 
perature T c , as compared to a noninteracting gas of the 
same density, as a function of an 1 / 3 , for small an 1 / 3 . 

In section II, the background is set by summarizing dif- 
ferent calculations and basic results. Section III contains 
the comments on various approaches to the T c problem. 
In the concluding section IV, I suggest some directions 
for additional calculations. 
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II. BACKGROUND AND PERSPECTIVE 

Early predictions for the shift in transition tempera- 
ture appear in the work of Lee, Yang and Huang, around 
forty years ago, as summarized in [1]. Recently, Huang 
has published [26] a prediction for an increase in T c due 
to the inter-particle repulsion, proportional to the square 
root of the interaction strength, T c — ~ V an 1 / 3 . 
Schakel in a short note [21], and Baym et al in [1], have 
both discussed the problem with this work; in particu- 
lar [1] contains an analysis showing how Huang's result 
is generated spuriously from a truncation of the virial 
expansion used at mean field level. 

A similar spurious result appearing in the literature 
is that of Toyoda [27], who in 1982 performed mean- 
field-like calculations deriving a Landau-Ginzburg type 
thermodynamic function. His prediction was a decrease 
of the transition temperature, AT C ~ — V an 1 / 3 . Again, 
this is believed to be an artifact of mean field approxi- 
mations. The present understanding is that there is no 
mean field effect on the transition temperature, and that 
higher-order correlations are required to find the leading 
shift in T c . 

A more recent work that also predicts a negative shift 
in T c is the canonical-ensemble work of Wilkens et al 
[20]. This work has been critiqued in [8]. I discuss some 
related issues in section III E. 

In addition, in view of the continuing appearance of 
negative-shift predictions, I give in section IIIB a simple 
argument for the increase of the transition temperature 
due to the addition of a repulsive interaction. This argu- 
ment makes use of eq (2), which is derived independent 
of any particular method of treating interaction effects. 

During the past several years, a large body of work 
has been published by Baym, Laloe and collaborators 
[1, 3-8], and a number of authors have done followup 
work [9-11]. These workers predict a positive shift in T c 
that to leading order is linear in the scattering length; 
AT C = cian 1 / 3 + 0(a 2 n 2/3 ). Attempts to calculate the 
coefficient ci, however, have continued to give fluctuat- 
ing results, and this has prompted the Baym group to 
predict [7] a logarithmic contribution at the next order, 
i.e., a —a 2 In a contribution. The presence of such a non- 
analyticity might explain the difficulty in any numerical 
estimate of c%. 
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Stoof's renormalization-group analysis [12, 13] of the 
three-dimensional bose gas has led to a prediction [13, 14] 
of an increase of the transition temperature proportional 
to a\ lna\, i.e., a non-analyticity at leading order. A lead- 
ing non-analytic behavior would actually explain even 
better the uncertain results obtained in attempts to cal- 
culate the coefficient c\ , based on the assumption of lin- 
ear shift AT C ~ cicm 1/3 . 

An older attack on the T c problem is that of Kanno 
[18, 19]. In a series of papers, Kanno introduced a "quasi- 
linear" canonical transformation approach to the many- 
body problem, which he used to calculate (among other 
things) the free energy and hence transition temperature 
of the bose gas. Since Kanno's work has mostly been ig- 
nored by the current T c community, I give a description 
of his work in section IIIF. I also show how some of his 
derivations can be performed in a more transparent way. 

Schakel's work on developing an effective theory [22] 
of the weakly non-ideal Bose gas has led to the predic- 
tion [23] of a shift in the transition temperature that is 
strictly linear in the interaction. The connection of this 
calculation with other approaches has not yet been ex- 
plored in detail. Some issues are pointed out in section 
HID. 

Our recent report [2] addresses the question of T c 
within a quasiparticle framework. Other than working 
out details of a quasiparticle description (ultraviolet di- 
vergence, effect of quasiparticle lifetime), we also provide 
a non-selfconsistent method of regularizing the infrared 
divergences appearing in the description of the transi- 
tion point. This results in a shift of T c approximately 
cx ay/\ loga|. The procedure is uncontrolled; however, 
this result points strongly to a non-analytic contribution 
in the dependence of T c on the interaction. 

Other relevant calcultaions include Kleincrt's "five- 
loop" determination [31] of the coefficient c\ mentioned 
above, assuming a puely linear leading shift, and the "lin- 
ear (5-expansion" calculations [32] of de Souza Cruz et al. 
Kleinert has questioned [31] the 5-expansion method, but 
a detailed analysis of the relationship of this approach to 
others seems to be lacking in the literature. 



A. Mean Field 

For a momentum-independent potential, the two mean 
field contributions to the self-energy (Hartree and Fock) 
are equal, S^^k, w) = Un + Un = 2Un. Since the self- 
energy is momentum- independent, it can simply be ab- 
sorbed into a redefinition of the chemical potential. 

« = Wo + ^9 3 / 2 (e' V ) 

Here // = fi — E^ 1 ) = ji — 2Un. The transition now 
simply corresponds to jj! = instead of /i = 0, and T c re- 



mains unchanged. The lesson here is that a momentum- 
dependent self-energy is required for a shift in T c . 

In second-order calculations, (e.g. in [2],) it is usual 
to have the lst-order self-energy already absorbed in the 
chemical potential, i.e., to use // as the chemical po- 
tential. This allows one to exclude those second-order 
diagrams which can be obtained by lst-order self-energy 
insertions into lst-order diagrams. 

Spatially extended interaction. In a real system, of 
course, the description of the interaction in terms of a 
single parameter (Uq or a) is valid only for very weak in- 
teractions or at very low temperatures/densities; if these 
conditions are not satisfied one has to take into account 
the detailed shape of the potential V(x) = Vqo-(x) which 
is not a delta function in real space (not a constant in 
momentum space) any more. In such a case we get a 
momentum-dependent Fock term, and hence a contribu- 
tion to T c from the first-order self-energy. Assuming a 
Gaussian V(x) of width a and height Vo, one gets a cor- 
rection to the spectrum that can be incorporated as an 
effective mass m* = m[l — Voa 2 nm/3]~ 1 , resulting in a 
negative shift [16]: 

AT C 1 2 

The negative shift at higher densities is present in liquid 
helium and presumably in other Bose fluids. However, 
for weak interactions, this shift is less important than 
the leading 0(a) effect, which is the current theoretical 
challenge. 

B. T c in Liquid Helium 

As a highly correlated liquid, the properties of liquid 
helium are not expected to be described well by the dilute 
weakly-interacting Bose gas model. However, the devia- 
tion of the lambda transition temperature (2.2K), from 
the BEC temperature of an ideal gas with the same den- 
sity (3. IK), can be mostly explained using a simple effec- 
tive mass description. In other words, the experimental 
T c is reproduced quite closely by replacing the helium 
atomic mass m by a quasiparticle effective mass m* in 
the ideal gas expression T c (0) = (27r/m[C(|)] 2/3 ) n 2/3 . 

The effective mass in the liquid exceeds the bare mass, 
because of the inertia of the medium which has to make 
way for any atom to move. This results in a lowering of 
T c . The effect is essentially equivalent to that described 
in the last section for a momentum-dependent interac- 
tion. 

Feynman [28] calculated a partition function for liq- 
uid helium that takes into account the effective mass en- 
hancement described above. Using this partition func- 
tion, Chapline calculated [29] a T c quite close to the ac- 
tual A-point of helium. This is equivalent to using the 
effective mass in the ideal gas expression for . 

According to current understanding, therefore, one ex- 
pects the following variation of T c with n^a: an initial 
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linear-like increase, followed by a maximum and then a 
decrease due to the effective-mass mechanism described 
above, until eventually AT C becomes negative. This kind 
of dependence has been seen in the Monte-Carlo calcu- 
lation of Griiter, Ceperley, Laloe [6], and in Reppy's ex- 
periments [30] described below. 

C. Relevant Experiments 

One could think of studying the T c problem experimen- 
tally, e.g., in the context of trapped atomic Bose gases. 
Unfortunately, trap effects dominate in the dependence 
of T c on the interaction for these systems. Turning on a 
repulsion for a trapped Bose gas causes an expansion of 
the gas, reducing the density and hence decreasing the 
transition temperature T c oc n 2//3 [33] . In contrast to the 
intrinsic effect under review here, trapping effects on T c 
are found at mean-field level already. 

A different experimental realization of the weakly non- 
ideal bose gas is the 4 He-Vycor system explored by Reppy 
et al [30]. Here the small-an 1 / 3 regime is reached by 
enormously decreasing the 4 He density. Reppy's investi- 
gation of the transition temperature shows a linear-like 
increase of T c at low an 1 / 3 , followed by a maximum and 
a decrease to values lower than Tc°^ . The general behav- 
ior is satisfying, but for the initial increase, the data is 
not y et g ood enough to choose between <~ a, <~ a In a, or 
<~ aVhi a-like behavior. 



III. COMMENTS ON VARIOUS APPROACHES 
A. Baym and collaborators 

The work of Baym, Laloe and collaborators [1, 3-8] es- 
tablishes definitely that a repulsive interaction increases 
T c , and that the dependence of AT C on the interaction is 
of linear order, as opposed to quadratic or <~ y/a. 

One should note, however, that the Baym group's work 
relies heavily on power-counting and dimensional argu- 
ments. This kind of argument typically does not catch 
logarithmic corrections. The possibility of leading non- 
analytic dependence of AT C on the interaction, therefore, 
needs to be further examined. 

It has been suggested [24] that the large-TV calcula- 
tion [4] is a reliable indicator that the leading shift is 
purely linear, because of intuition gained from other sys- 
tems where large- TV calculations capture all logarithms. I 
would like to point out, however, that the TV ^ 2 models 
belong to a different universality class, and have differ- 
ent excitation spectra as compared to the usual (TV = 2) 
Bose gas. More care is therefore required in interpreting 
results from the large- TV approach. 

The Baym-Laloe collaboration seems to have been the 
first to relate the single-particle spectrum at T c to the 
shift AT C . Another noteworthy aspect is the use of Ursell 
operators [1, 5], and a corresponding diagram technique, 



to describe particle correlations. (As far as I can tell, 
for this problem, the Ursell operator technique does not 
give any extra information compared to the usual Green 
function techniques.) 

B. Direction of Shift, Infrared Spectrum 

Once the shift in T c is related to the modification of 
the single-particle spectrum, the direction of the shift can 
be argued using well-known facts about the spectrum of 
the Bose gas. 

We will use the following expression for the shift, as 
appearing in [2]: 

Here / c is the Bose function at the critical point; f c (x) = 
[e x / kTc — and £k is the quasiparticle spectrum at 

Tc i.e., £ k = e k + E(k,& ; T = T c ). 

The expression used by Baym et al [1, 3], AT C cx 
J dkU(k)/[k 2 + U(k)}, is the high-temperature (zero Mat- 
subara frequency) version of (2). Noting that the major 
contribution to (2) comes from the infrared, one can write 
approximately 

/•cutoff 

AT c cx / dkk 2 ^ 1 -^ 1 ] . (3) 
Jo 

This equation might introduce or leave out logarith- 
mic corrections; however the following argument is un- 
changed for (2). 

If the quasiparticle spectrum £k is "harder" than the 
bare spectrum e k = k 2 /2m, i.e., if £k is sub-quadratic 
in k or quadratic with an effective mass to* < m, then 
the ek integral dominates and the shift is positive. On 
the other hand if the spectrum at the transition point is 
"softened" by the interaction, then eq (3) [ or (2)] leads 
to a decrease in the transition temperature. 

Since the weakly interacting bose gas is known to have 
a linear (Bogoliubov) infrared spectrum below the critical 
temperature, and a quadratic spectrum above T c , it is 
extremely reasonable to assume, even without any RG 
or perturbation-theoretic calculations, that the spectrum 
at the critical point should be something intermediate 
between linear and quadratic. As pointed out above, this 
means a positive shift in the transition temperature. 

The long-wavelength spectrum at the transition point 
should be £ k ~ k 2 ~ n , where r\ is the anomalous dimen- 
sion. The present system falls in the same universal- 
ity class as the XY model, or the TV = 2 quantum ro- 
tor model. For this universality class, large- TV and e- 
cxpansion calculations give respectively rj ss 0.14 and 
77 w 0.02, and a recent Monte Carlo calculation [35] has 
produced the value 77 w 0.038. 

Second-order treatments of the self-energy (e.g., the 
Baym group's [1, 3] and our [2] work) over-estimates the 
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modification of the spectrum, giving a k 3 / 2 infrared be- 
havior (rj = 0.5). On the other hand, some other treat- 
ments (such as Stoof's [13] and Schakel's [22, 23]) neglect 
the modification altogether (77 = 0). 



C. Renormalization group treatment (Stoof) 

The RG flow equations derived by Bijlsma & Stoof 
[13] allow one to numerically locate the transition point. 
Under RG flow, the physical effective chemical potential 
flows to positive and negative values, if one starts out 
respectively in the Bose-condensed and the symmetry- 
unbroken phase. The critical point is therefore identified 
by finding parameters for which the chemical potential 
flows to zero. Stoof finds [13, 14] that the dependence 
of AT C on a I A has more structure than purely linear; 
AT C oc (a/A) ln(a/A). 

The differences between Stoof's and the Baym group's 
results merit further study. Stoof has suggested [14] that 
the source of the difference is not the infrared but the 
ultraviolet. The correction to the spectrum, U(k) oc 
S(k, Ek) — S(0, 0), has an ultraviolet In k behavior in the 
Baym group's treatment [1, 3]. In the RG framework, 
the corresponding quantity is the running coupling con- 
stant /x(A), which shows no ultraviolet logarithms in the 
calcultation of ref [13]. This leads to a logarithmic dif- 
ference in the two predictions for T c . Stoof also suggests 
[14] that an improved RG calculation, including the mo- 
mentum dependence of the effective two-body interaction 
at the transition point, would shed light on the difference 
between the two approaches. 



D. Schakel: IR regularization 

In treating his effective-action theory near the transi- 
tion point, Schakel [22, 23] is also faced with infrared (IR) 
divergences. One of the remarkable features of this work 
is an innovative analytic continuation procedure used to 
regulate IR problems. 

Baym et al have remarked [1] that Schakel's prescrip- 
tion ((-function regularization) must be incorrect be- 
cause it gives an incorrect result for the compressibility 
of the ideal gas. This argument, unfortunately, is derived 
from an equation that Schakel derived explicitly for the 
nontrivial broken-symmetry ground state [23, 24] of an 
interacting gas (fi > 0). The same equations are not ap- 
plicable to the noninteracting case (fi < 0) for which the 
structure of the vacuum is different (more trivial) . The 
situation is therefore more complicated than indicated in 
Baym et aFs note. In response, Schakel has provided a 
more detailed explanation and justification [25] for the 
(■-function continuation procedure. 

It seems to me that the regularization procedure it- 
self is quite possibly sound. My concern with Schakel's 
derivation is different — just below the transition point 



(broken-symmetry phase) , Schakel uses a quadratic spec- 
trum. A more accurate treatment would use the phonon- 
like spectrum E^ s=a yj e 2 . + 2Un a (T)e] i . The effect such 
a spectrum would have on the thermodynamic potential 
or on T c is not clear. (Redoing Schakel's T c calculation 
with a phonon-like spectrum seems quite difficult.) 



E. Canonical Ensemble: Role of No(T c ). 

In Wilkcns et aFs paper [20], calculations have been 
performed using a canonical ensemble of N particles, and 
the shift AT C obtained is linear and negative in the scat- 
tering length. It has been pointed out [8] that the er- 
ror lies in evaluating quantities (derivatives) at No = 0, 
instead of evaluating them at the expectation value of 
No(T c ). Mueller et al [8] have also used finite-size scaling 
ideas in combination with a first-order calculation in the 
canonical ensemble, and obtain a positive linear shift. 

The idea that the zero-momentum state occupancy Nq 
becomes "microscopic" at and above the transition needs 
to be clarified. One finds lim N ^ (N /N) = for T > 
T c , where N is the total number of particles, but not 
necessarily ~ 0{l). In fact, for both the ideal and 
the interacting Bose gas, N (T C ) <~ N 2 / 3 , which is large 
compared to 1 though small compared to N. 

In view of the continuing appearance of the No(T c ) ~ 
O(l) misconception, I give here a simple derivation of 
No(T c ) = ^yN 2 ^ 3 , together with an expression for the 
coefficient 7, for the ideal-gas case. 

Let us consider the ideal-gas relationship N — N + 
^■93/2 (e a ), where a — —(3fi, and g v refers to a Bosc- 
Einstein integral function [37]. Since a is small near the 
transition, one can use the Robinson expansion [37] for 
g 3 / 2 - Also, inverting No = l/(e /3£o+Q — 1) one gets a w 
1/Nq. Combining these yields, near the transition, 



iv 3 / 2 + 



V 



.V„-^(2V^) 







Since nA 3 = ((|) at the transition, the coefficient of \/Nq 



needs to vanish. Therefore 



N (T C ) 



am 



2/3 



(4) 



The k = occupancy No can actually be calculated for 
all T, using an extra Lagrange multiplier /io conjugate 
to No , in addition to usual chemical potential \i coupling 
to N. Using this "two chemical potentials" formalism 
[34], one can show that for the ideal gas (No) — 1 occurs 
at a temperature T w 2.6 Tc°\ well above the transition 
point. 

For the interacting-gas case, an equation correspond- 
ing to (4) is difficult to derive, but finite-size scaling ar- 
guments [8, 17] show that N ~ N 2 / 3 continues to hold, 
even though the nonideal Bose gas belongs to a different 
universality class. 
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The definition of the transition point used in ref [20] is 
in terms of the probability 



P n (T,N) cxtr 



6 - e~ 0H 8 - 

u N ,n c N,N 



(5) 



for having a condensate occupancy No = n. The claim 
is that P n (T, N) decreases monotonically with n for all 
temperatures above the transition, and has a peak at 
nonzero n for in the condensed phase, and therefore the 
transition should be defined as the thermodynamic limit 
of the zero of D(T, N) = P Q (T, N) - P^T, N). 

The correspondence between the canonical particle- 
counting statistics, and the grand-canonical ideas of av- 
erage occupancies, is not trivial. One important point to 
note is that a monotonically decreasing P n does not im- 
ply a zero expectation value for No, even though P n has 
a maximum at n = 0. (For example, using P n oc e~ an 
gives {No} — J2 n n ^ n = V[ eQ — !]■) The stipulation of a 
monotonically decreasing P n therefore does not imply a 
zero or 0(1) value of N at the transition. 



F. Kanno's Approach: Using Modified 
Thermodynamic Potential 

Kanno developed an unusual approach [18] to study 
many-body systems, and also used this approach to cal- 
culate Bose gas properties [19]. 

The formalism. 
Kanno's calculation of thermodynamic quantities in- 
volves a "quasilinear" canonical transformation [18] of 
the free-particle operators b {b'}: 



to point out that many of his physical results, if not all, 
are equivalent to results obtained in more mainstream 
approaches. For example, the quantities A{p) and B{p), 
appearing in the second paper of the series [18], are equiv- 
alent to the normal and anomalous self-energies that ap- 
pear in modern descriptions. 

The shift in T c . 
Kanno's defines [19] a modified thermodynamic potential 

O'(Ar ,M',T) = -ilntr'e-^-^') , 

with an unusual Lagrange multiplier jj conjugate to the 
number of non-condensate particles, TV' = N — No, as op- 
posed to the usual chemical potential which is conjugate 
to the total particle number N. It follows that 



N-No = N' 



dn'(No,n',T) 



(7) 



Once Q' is calculated in some approximation, comparing 
this relation to N — N n = V gzii{ e } / 'A 3 for the free bose 
gas at the transition point yields the shift AT C , in terms 
of the quantity /i'(T c ). 

The quantity // itself is obtained from a consideration 
of the free energy F(N ,N,T} = Q'{Nq, /j,', T) + fi'N'. 
The actual value of iV chosen by the system is the one 
that minimizes F{Nq, N,T), i.e., 



= 



dF 



dNo 



dSl' 



T=T C 



8No 



(8) 



T=T C 



k = it + ^2 k, l}i]ikii 

j;k,l 

+ Ai ti> h > l > m > n )l]lllllmln + ■■■ , (6) 

resulting in a Hamiltonian having a particular (diagonal- 
like) form. The formalism involves examining the proper- 
ties of the coefficients A, and the coefficients appearing in 
the quasi-diagonal hamiltonian. Applying this formalism 
to the weakly repulsive bose gas [19], one can derive the 
Bogoliubov spectrum, entropy and pressure above and 
below the transition, and many other quantities. 

Kanno does not provide direct connections to the usual 
language of perturbative many-body theory. I would like 



This equation allows one to express ^'(T c ) in terms of 
the interaction strength. 

Computing O' in more usual language. 
As advertised in section II, I show here how the modi- 
fied thermodynamic potential f2' (which is the quantity 
relevant to the T c problem) can be obtained from a dia- 
grammatic second-order perturbative calculation. 

The Feynman rules for a diagrammatic perturba- 
tive calculation of the modified therodynamic poten- 
tial Q'{fj,',No,T} are very similar to those for the usual 
£l(p,T). One has to use a momentum-dependent chem- 
ical potential /ik = //(l — <5k,o)) and separate out k = 
from the momentum-sums at each stage. Two first order 
and two second-order diagrams need to be calculated: 




Other second-order diagrams are excluded by considering the chemical potential to be renormalizcd upto first order, 
as discussed in section II A. 
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The zeroth and first order terms can be evaluated exactly. UV divergences appearing at second order can be 
removed by eliminating the bare potential U in favor of the 2-body t-marix, as done for the self-energy in [2]. The 
result is 



V / p / q e P -eq-e P -q 



+ 2t 2y I I I /(e p -M0/(6 k -M0/(^-^) (Q) 

Jp Jq Jk e p + £ k — e q — e p+k-q 



The expression for O' obtained in Ref [19] is essentially 
the same as eq (9). 

Comparison & remarks. 
The parameter fJ,'(T c ) corresponds to the infrared cut- 
off used in the the non-selfconsistent regulariztion of [2] . 
Eq (8) gives the same result for the parameter n'(T c ) as 
that obtained in [2] for the chemical potential \x using 
\x = £(0,0). The reason for this is that dF/dN repre- 
sents the energy of a zero-momentum quasiparticle with 
respect to the chemical potential, i.e., E(0, 0) — fj,. There- 
fore eq (8) is equivalent to the usual condition S(0, 0) = yU 
for the Bose-condensation point. 

The shift AT C obtained using eqs (7), (8), (9) is the 
same as that obtained by our non-selfconsistent regular- 
ization scheme in [2] . 

The distinctive idea in Kanno's work is the introduc- 
tion of a modified thermodynamic function f2' such that 
—dCl'/d/j,' replaces the right side in the criticality con- 
dition N = (V7A 3 )C(§). This bypasses the issue of self- 
consistency and the modification of the spectrum at T c , 
and the use of the chemical potential as a cutoff folllows 
naturally. 



IV. SUMMARY, FURTHER CALCULATIONS 

To summarise, this article reviews aspects of the calcu- 
lation of the T c shift due to a small repulsive interaction 
in a Bose gas. This question has generated a number 
of calculations involving diverse techniques; the present 
discussion is a first step toward developing an organized 
overview of the different approaches. 

The issue that needs to be resolved urgently, in my 
opinion, is the possibility of a logarithm-like correction 
to the leading shift in T c . Presumably, a better under- 
standing of the spectrum at the critical point (both at 
the infrared and the ultraviolet) should help settle this 
issue. Two directions need to pursued in this regard. 

First, as suggested by Stoof, an improved RG calcu- 
lation, taking into account the momentum dependence 
of the effective interaction, would improve the ultravio- 
let behavior of the spectrum in his formalism. Also, it 
might be worth retaining the wavefunction renormaliza- 
tion (critical index 77), which changes the infrared spec- 
trum at T c from ~ k 2 to ~ k 2 ~ v . This is neglected in 
[13] on the grounds that 77 is small. 



A second approach suggests itself when we note that 
the quasiparticle treatments [1-3] have all calculated self- 
energies in terms of the two-body t- matrix t = Aira/m. 
This quantity describes all possible collisions between two 
particles, without taking into account the fact that the 
surrounding medium has an effect on these collisions. 
Many-body corrections arising from the surrounding gas 
can be treated by using the many-body T-matrix, which 
is the sum of ladder diagrams. The two-body t-matrix 
has no temperature-dependence, while the temperature- 
dependence of the effective potential (vanishing at T c ) is 
captured by the many-body T-matrix [13, 15, 36]. Stoof 
suggests [14] that this might explain why the quasiparti- 
cle analysis (perturbative treatment of spectrum modifi- 
cation) over-adjusts the infrared spectrum. 

Determining the self-energies and quasiparticle spec- 
tra in terms of the many-body T-matrix has been at- 
tempted [34] , using the relations provided by Shi & Grif- 
fin [36] between the bare potential U, the t-matrix, and 
the many-body T-matrix. Our attempt was unsuccessful 
due to the difficulty of disentangling multiple frequency 
and momentum inter-dependences at second-order. 

Finally, I mention some other work that would be in- 
teresting to the T c community. 

• One possibility is to use finite-size scaling ideas, in 
conjunction with second-order perturbation theory 
in the grand-canonical ensemble. This would be 
a complementary to the calculation of Mueller & 
Baym [8], which uses finite-size scaling with first- 
order perturbation theory in the canonical ensem- 
ble. 

• It would be nice to see Schakel's calculation of the 
thermodynamic potential, just below T c , done with 
the spectrum treated more properly, either with 
a Bogoliubov-like spectrum or with the E^(T C ) <~ 
k 2 ~ v spectrum. 

• As mentioned in section II, a detailed comparison 
is required between the 5-expansion work [32] and 
other approaches to T c . 
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